Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences

The spatial resolution of optical imaging systems is fundamentally restricted by the diffraction limit. However, in widefield live-cell microscopy, the achievable resolution is further constrained by the specimen motion, which indicates the existence of a fundamental spatio-temporal resolution trade-off between signal accumulation during the full frame integration and the resulting motion blur. To improve the fidelity with which moving objects can be imaged, a quantitative understanding of this spatio-temporal trade-off is necessary. Here, we present a systematic analysis of motion-induced resolution dynamics measured with spectral signal-to-noise ratio (SSNR). We developed a simulation framework which models the image formation of objects undergoing arbitrary motion, to evaluate the degradation of the spatial resolution under translational and rotational dynamics. Our results demonstrate that for translating objects, the spatial resolution is anisotropically reduced as a function of the orientation of the object relative to the motion vector, leading to the spectral signal-to-noise ratio degrading by up to 50% and the resolution by up to 40% for a 90{degrees} change in the motion direction. Furthermore, we show that for rotational motion, conventional radially averaged metrics such as the Fourier Ring Correlation are not able to quantify the effects of angular blur. On the other hand, the SSNR is able to accurately quantify this degradation. These findings underscore the necessity of an object-oriented imaging approach, in which acquisition parameters such as exposure time are tuned to specific biological spatio-temporal characteristics to optimize the trade-off between motion blur and spatial fidelity.

Uncertainty quantification (UQ) in computational heart models is essential for reliable cardiac digital twins (DTs) in personalized medicine, yet remains challenging. Traditional Monte Carlo and stochastic Galerkin methods often become impractical in the high-dimensional, nonlinear state variable and parameter spaces of cardiac electrophysiology and mechanics. This article introduces a framework for learning a joint probability density over cardiac observables and model parameters, enabling the characterization of statistical dependencies across a large number of variables in patient-specific cardiac DTs. By sampling from this density and conditioning on available data, useful predictive distributions can be constructed, allowing uncertainty to be propagated through the model and quantified in terms of variability. Conditional regression can then be performed directly on this learned density, enabling systematic exploration of interdependencies among observables for both predictive inference and model design. The statistical methodology adopts a geometry-aware generative learning framework, recently introduced by the authors, that decouples the learning of data geometry from sampling. First it identifies a low-dimensional latent representation that captures the intrinsic structure of the data and its multiscale geometric features. A stochastic differential equation is then formulated directly in the low-dimensional latent space to generate samples efficiently; these are subsequently mapped back to the high-dimensional space of cardiac states and parameters through a smooth lifting operator. We demonstrate the approach on a ventricular arrhythmia prediction benchmark, where the learned joint probability density enables the construction of predictive distributions over key parameters (e.g., conductivities, fibrosis patterns) through sampling and conditioning. This enables uncertainty to be propagated and quantified through sampling and conditioning on the learned joint density, with substantially fewer model evaluations than conventional UQ methods.

We introduce a workflow to identify oligomeric structures that are recorded with single-molecule localization microscopy (SMLM) under cryogenic conditions. Typically, these oligomers are assumed to consist of protomers arranged as equilateral two-dimensional polygons and every protomer is labeled with a dye molecule for visualization. Unlike previous work, we consider scenarios in which the sample plane has an unknown orientation relative to the focal plane. Our contribution is a high-precision plane-fitting algorithm to determine the sample plane, combined with geometrical transformations and two circle-fitting algorithms to identify the oligomeric structures. Our simulations on synthetic data demonstrate that the proposed workflow achieves high accuracy in estimating both the unknown tilted plane and the oligomer size.

Quantum kernel methods offer a potential advantage for classification tasks in high-dimensional feature spaces, yet their practical benefit critically depends on how input features are prepared. We compare five dimensionality reduction strategies--principal component analysis (PCA), Gaussian random projection (RP Gaussian), sparse random projection (RP Sparse), partial least squares (PLS), and uniform manifold approximation and projection (UMAP) -- as pre-processing steps for quantum kernel support vector machines (SVMs) applied to trabecular bone classification from synthetic micro-computed tomography (micro-CT) data. Using a custom procedural generator based on Gaussian random field zero-crossings, we produced 500 synthetic trabecular bone volumes with controlled morphometric properties such as bone volume fraction (BV/TV), trabecular thickness (Tb.Th), number (Tb.N) and spacing (Tb.Sp). Texture features extracted from grayscale slices are reduced to 8-dimensional quantum circuit inputs via each method, then classified using both classical radial basis function (RBF)-SVMs and quantum kernel SVMs with ZZ feature maps on a statevector simulator, both evaluated with 5 x 5 repeated stratified cross-validation (25 folds). Our results show that UMAP is the only reduction method where the quantum kernel remains competitive with the classical baseline. Under repeated cross-validation, UMAP showed a +0.032 accuracy gap favouring the quantum kernel (Dietterich 5 x 2 CV p = 0.177); however, validation on 10 fully independent datasets--each with independently generated samples, separate reduction fits, and separate kernel matrices -- reversed the sign to -0.030 (paired t-test p = 0.123; Wilcoxon p = 0.193; quantum wins 3/10 datasets), indicating that the apparent advantage was likely inflated by fold dependence. Nevertheless, UMAPs gap remains small and non-significant in both analyses, whereas all linear methods (PCA, RP Gaussian, PLS) show substantial quantum deficits of -0.090 to -0.116 across BV/TV classification, with PCA and PLS remaining significant under corrected tests (5 x 2 CV p = 0.004 and p = 0.007 respectively). We additionally evaluate quantum kernel ridge regression for continuous morphometric prediction, finding that ZZ quantum kernels fail uniformly at regression (negative R2 for all methods except PLS at 4 qubits), suggesting that the ZZ kernel captures decision boundaries but not smooth metric structure. These findings provide practical guidance for feature engineering in near-term quantum machine learning pipelines and demonstrate that the choice of dimensionality reduction can determine whether quantum kernels remain competitive with classical baselines.

The Poisson-Nernst-Planck (PNP) system is an accurate model of electrodiffusion of ionic species. It is commonly used in situations where nanoscale resolution is required, for instance close to ion channels in the membranes of biological cells. The inherent stiffness of the equations has made them challenging to solve and has limited the applicability of the system. In particular, the time step required for stable solutions has typically needed to be very short (nanoseconds), which makes simulations on the time scale of an action potential (milliseconds) difficult. Recently, it has been observed that avoiding operator splitting and instead solving the concentration equations and the electrostatic equation in a coupled manner relaxes the time-step limitation considerably. However, no theoretical explanation of this observation has been provided. Here, we aim to explain why the coupled scheme allows much larger time steps. We illustrate the mechanism by considering special cases that define necessary, but not sufficient, conditions for stability. We also show that these conditions remain relevant for the fully coupled PNP model in 3D.

A major challenge in neuroscience is to predict how currents in nanodomains affect voltage and ionic concentrations. Cable and Rall theory provide analytic current-voltage relations by neglecting concentration gradients, and the impact of concentration gradients is usually studied numerically with the Poisson-Nernst-Planck (PNP) model. A precise quantitative understanding of the combined dynamics remains limited because analytic current-voltage-concentration relations are missing. In this work we derive such relations using a novel approach based on cross-diffusion equations. For narrow cylindrical domains, we derive time-dependent and steady-state expressions that explicitly show how currents affect voltage and ionic concentrations. We find that the influx of only one ion can significantly change the concentrations of all the other ions even if no channels for these ions are present. After a current injection we compute a biphasic voltage transient where the small-time asymptotic corresponds to the steady-state solution of the cable equation. We show that the accuracy of cable theory prediction for the voltage depends on how the current is distributed among the various ions. Finally, we develop an iterative method to accurately compute steady-state profiles for voltage and concentrations using first-order results by subdividing a cylinder into small segments.

Quantifying the diagnostic dispersion of inferred parameter distributions is a challenge in uncertainty-aware modeling. Scalar summaries such as credible interval width are topology-blind; fundamentally different posterior morphologies can yield identical scores, obscuring whether a parameter is precisely estimated or constrained to a range. We propose a Composite Certainty Framework that addresses this metric degeneracy by aggregating five complementary uncertainty metrics including interquartile range, standard deviation, full width at half maximum, Shannon entropy, and mass width. These metrics are aggregated through non-parametric Borda rank voting into a single, unitless consensus certainty score. Applied to a simulation-based inference pipeline for a finite-element model of the human middle ear tuned to cadaveric acoustic measurements, the framework reveals parameter-specific identifiability profiles invisible to any individual metric. It produces two actionable clinical thresholds: (1) the maximum tolerable measurement noise for reliable parameter recovery, and (2) the minimum simulation budget for posterior convergence. We demonstrated that no single metric captures all aspects of posterior dispersion, as spread-based metrics and entropy diverge systematically for clinically critical parameters, whereas their aggregation produces a consensus reflecting genuine diagnostic certainty. The framework is generalizable to any model-based diagnostic pipeline where posterior distribution not merely its coverage, but determines clinical certainty.

For decades, the photon counting histogram (PCH) was used as the sole method to quantify fluorophore numbers in a diffraction-limited focal volume. This technique combines spatial excitation profiles, and the distribution of photon counts to register the photon emission statistics of individual fluorophores. However, this approach has not yet been transferred to widefield fluorescent imaging due to the lack of fast and single photon sensitive camera sensors which can capture the photon emission statistics of a single fluorophore. Here, we explore avenues towards quantitative analysis of the active fluorophore number by leveraging recent advancements in single photon avalanche diode (SPAD) array technology. Binary exposures of a SPAD array can be synchronized with picosecond laser pulses to measure the PCH in a widefield setting. Then, by modeling the statistical relationship between the active fluorophore number and the PCH in a region of interest following a laser pulse, we can perform Bayesian inference of this number. The model is demonstrated experimentally by counting quantum dots and various numbers of fluorescent dye molecules bound to DNA origamis. We find that this method has several important applications in widefield microscopy, including enhanced localization microscopy and constrained fitting of multiple unresolvable fluorescent emitters.

This study presents a proof-of-concept cyber-physical architecture integrating a SEIR epidemiological model (Susceptible-Exposed-Infectious-Recovered), implemented in MATLAB, with a simulated Internet of Things (IoT) acquisition and transmission stage based on the ESP32 microcontroller and the ThingSpeak platform. The system generates synthetic biomedical signals of body temperature and peripheral oxygen saturation (SpO2), structured across three levels: circadian variation, scheduled pathological episodes, and Gaussian noise. These signals feed a dual parametric coupling function that dynamically updates the SEIR transmission parameter as a combined function of body temperature and oxygen saturation deviations from their clinical reference values. The proposed architecture is organized into four functional phases: measurement, communication, computational processing, and feedback. Five simulated clinical scenarios were evaluated, ranging from normal conditions (T = 36.5 {degrees}C, SpO2 = 97%) to fever with severe hypoxia (T = 38.5 {degrees}C, SpO2 = 88%), yielding basic reproduction number (R0) values between 4.20 and 5.38, and peak infected proportions between 29.9% and 35.2% of the simulated population (N = 1,000). A sensitivity analysis on the coupling coefficients, with {+/-}50% variation from nominal values, showed that the oxygen saturation coefficient is the most influential parameter on R0 (range = 0.76) compared to the thermal coefficient (range = 0.42), with monotonic and predictable behavior across the entire evaluated parametric space. The primary contribution of this work is system integration: we propose a reproducible platform connecting biomedical simulation, IoT communication, and epidemiological modeling through parametric coupling in a controlled environment. All data used are entirely synthetic; a retrospective calibration with real Colombian data from the first epidemic wave of 2020 confirmed the epidemiological consistency of the model, with a calibrated R0 of 1.85 and a Pearson correlation of 0.930. Results should be interpreted as evidence of architectural feasibility, not as clinical or epidemiological validation. Author SummaryThe COVID-19 pandemic made it clear that epidemiological surveillance systems need tools that combine accessible technology with mathematical models capable of anticipating disease spread. In this work, we built a proof-of-concept platform connecting three elements: a low-cost electronic sensor based on the ESP32 microcontroller, a cloud communication platform (ThingSpeak), and a mathematical model that simulates how an epidemic spreads through a population. The sensor generates synthetic data on body temperature and oxygen saturation that, through a mathematical formula we designed, dynamically modify the rate of contagion in the model. We evaluated five clinical scenarios, ranging from normal conditions to fever with severe hypoxia, and analyzed how sensitive the results are to changes in the system parameters. We found that oxygen saturation has a greater influence on the estimated contagion potential than body temperature. Although all data are synthetic, this platform demonstrates that it is possible to integrate low-cost sensors with epidemiological models in real time, opening a viable pathway for early warning systems in resource-limited settings.

The purpose of electronic disease syndromic surveillance (EDSyS) systems is to detect hazardous pathogens and other unusual signals in health surveillance data before such events are identified by an individual clinician or healthcare facility. However, EDSyS systems have primarily been evaluated using simulated health surveillance data, which do not necessarily capture the richness and complexities of real-world health data. We have updated and extended an existing EDSyS system, EpiDefend, which combines ensemble forecasting and recursive Bayesian estimation in a particle filter framework that supports demographic and spatial structure. We simulated the release of several pathogens, both infectious and non-infectious, and injected the resulting cases into a real-world health data set. Here we evaluate EpiDefend's sensitivity and specificity in detecting these simulated releases, and measure the time to detection against pathogen-specific estimates of the time to clinical detection, as informed by clinicians and microbiologists. We show that for diseases where clinical diagnosis can be challenging, such as Q fever (Coxiella burnetii) and tularaemia (Francisella tularensis), EpiDefend can reliably beat the time to clinical detection. In contrast, for pathogens that can be clinically diagnosed relatively quickly, such as inhalational anthrax and pneumonic plague, it is extremely difficult to beat the time to clinical detection. Our results suggest that EpiDefend may be able to reliably detect real-world introductions or releases of some pathogens at low false-alarm rates before a clinical diagnosis would be confirmed, and this would represent a landmark achievement for EDSyS systems.

We systematically evaluate and compare the performance of classical statistical methods (ARIMA), mechanistic compartmental models (SEIR), modern deep learning architectures (LSTM, DLinear, Autoformer), and an emerging time-series foundation model (TabPFN-TS) to forecasts the incidence of Influenza-Like Illness (ILI) across nine European countries. The models are benchmarked against a naive baseline and a multi-model ensemble (RespiCast) created by an initiative of the ECDC. In line with the operational practice of existing forecasting hubs, our entire evaluation is explicitly optimized for short-term horizons (1 to 4 weeks ahead). Interestingly, we found that the foundation model TabPFN-TS allows for great zero-shot inference capabilities. Without any task-specific retraining, it successfully overcomes extreme data scarcity to consistently outperform all other individual architectures, frequently rivalling or surpassing the RespiCast ensemble. Our results highlight how deep learning architectures are severely constrained by extreme data scarcity, typical in epidemic forecasting, requiring targeted endogenous data augmentation to reduce predictive errors. Within the deep learning class of models, we observe that simpler architectures (such as DLinear and LSTM) frequently exhibit greater robustness and outperform complex, attention-based models (such as Autoformer) when data is constrained. Finally, our results show how a weighted ensemble, constructed by fusing all the models, delivers highly robust forecasts in all regions considered. Overall, our findings showcase the transformative potential of zero-shot foundation models in epidemic forecasting and confirm the importance of multi-model ensembles.

Phylodynamics bridges the gap between epidemiology and pathogen genetic data by estimating epidemiological parameters from time-scaled pathogen phylogenies. Multi-type birth-death (MTBD) models are phylodynamic analogies of compartmental models in classical epidemiology. They serve to infer the average number of secondary infections R and the infection duration d. Moreover, more complex MTBD models add extra parameters, such as the average length of the incubation period or the proportion of superspreaders in the infected population. However, these additional parameters come at an important computational cost: Apart from the simplest, BD, model, MTBD models do not have a closed-form solution and require numerical methods for their likelihood computation. This leads to increased computational times and potential numerical errors. Therefore, the BD model remains the favorite researchers choice for real dataset analyses, and is often applied even in cases where more complex epidemiological aspects are present. We investigated, using simulations, how model misspecification influences inference of R and d in the phylodynamic framework. We showed that the use of models not accounting for various epidemiological aspects leads to bias. In particular the simplest, BD, estimator tends to underestimate R in the presence of super-spreading or incubation, which might be dangerous from the public health prospective. However, deep-learning-based estimators for complex models, which account for multiple epidemiological factors, perform well both on the data where those factors are present and where they are absent. This advocates for the use of complex epidemiologically realistic estimators, whose design has recently become possible thanks to deep learning.

Infectious disease time series often show signs of epidemic transitions, such as the peaks and troughs of the time series. In these time series, key system parameters can lead to catastrophic changes in the dynamical system behaviour (often called critical transitions). Modellers have increasingly shown that early warning signals can anticipate these transitions, both critical and non-critical, in infectious disease time series. Existing methods, however, generally focus on univariate time series data, or ignore spatiotemporal patterns that may be present as a disease spreads through a population. Recent ecological literature developments expand existing temporal and spatial methods to consider the covariance matrix of multiple, related time series. However, many of these proposed signals still make an assumption of stationary time series/system equilibrium. Whilst often true in ecological modelling, disease systems are seldom at equilibrium. In this paper, we propose the usage of the eigendecomposition of the non-stationary covariance matrix as a more suitable early warning signal for epidemiological data. We first analyse the expected trends in the eigenvalues and eigenbasis of the covariance matrix on approach to a transition. Next we apply these methods to a spatially-structured susceptible-infectious-recovered model to explore how the eigenbasis may provide extra information to modellers. Finally, we test these methods on SARS-CoV-2 case data during the 2020-2021 pandemic period in England.

Understanding how neuronal populations interact to encode and transform sensory information is a fundamental challenge in computational neuroscience. Most existing studies, however, study neural encoding, behavioral readout, and functional connectivity as disjoint problems. Two-photon calcium imaging enables simultaneous recording of large neuronal ensembles in vivo, driven by diverse stimuli and eliciting distinct behaviors. However, extracting directional functional connectivity metrics as well as encoding and readout properties of neurons from such data remains difficult due to indirect and noisy observations of spiking activity, slow temporal dynamics, and the latent interplay between external stimuli and endogenous neural processes. Here, we introduce a unified conceptual and operational modeling and inference framework for directly extracting functional Granger causal (GC) effects between neurons, from external stimuli to neurons, and from neurons to behavior, from two-photon imaging data, in the sense of Granger. Inspired by the intersection information framework, we also identify neurons that encode features of sensory stimuli that inform behavioral readout. The resulting GC networks together with the taxonomy of functional sensori-behavioral relevance, which we call G-taxonomy, provides a powerful statistical analysis framework, enabled by the integration of several techniques including state-space modeling and inference, variational inference, and point processes. We applied the proposed framework to simulated and experimentally-recorded two-photon imaging from the mouse auditory cortex (A1) during both passive listening and active tone discrimination. Our simulation studies reveal significant improvement of our proposed methodology over existing techniques. Analysis of experimental data from the mouse A1 identifies distinct groups of cells with diverse sensori-behavioral relevance, as well as changes in functional connectivity associated with correct vs. incorrect behavior. In summary, this work provides a principled and data-driven methodology for uncovering directional interactions among the neurons, sensory stimuli, and behavior, all within the same statistical framework, offering new insights into how distributed cortical populations transform sensory inputs into behaviorally relevant representations. Author SummaryThe brain processes sensory inputs through the coordinated activity of large networks of neurons and produces readouts that elicit behavior. Understanding how information flows and is processed through these networks is a central goal of neuroscience. In this study, we present a new computational framework that identifies directional interactions among neurons in an ensemble as well as from sensory stimuli to neurons and from neurons to behavior. Utilizing the Granger formalism to identify directional effects, as opposed to common correlational measures, our framework extracts said effects directly from two-photon calcium imaging data. We tested our proposed method on both simulated data and recordings from the auditory cortex of mice during passive listening and active tone discrimination tasks. Our method revealed diverse groups of neurons in the auditory cortex with distinct functional roles and relevance to sensori-behavioral integration. Our framework provides a new way to study the flow of information in the brain and can be broadly applied to uncover neural computations across sensory and cognitive systems.

Following ER Ca2+ depletion, Ca2+ release-activated Ca2+ (CRAC) channels are activated by STIM1 at ER-plasma membrane junctions. The restricted localization and low conductance of the CRAC channel (<40 fS) precludes single-channel recordings, limiting studies of CRAC channel gating. Here we describe an optical approach to characterize the gating of HaloTag-fused Orai1 channels labeled with JF646-BAPTA, a Ca2+-sensitive fluorescent dye. While Ca2+ influx through single channels generates fluorescence fluctuations, identifying true gating events is complicated by stochastic transitions of JF646-BAPTA to a non-fluorescent state. To overcome this, we combine TIRF microscopy with whole-cell voltage clamp to control the driving force for Ca2+ entry. We show the open channel intensity at -100 mV reflects Ca2+ saturation of the dyes on each channel, while the closed-channel intensity is defined by the fluorescence at +30 mV, where influx is absent. True gating events can be identified from transitions between the open- and closed-channel levels, distinguishing them from transitions to a non-fluorescent state. We describe the gating behavior of CRAC channels activated by STIM1 after store depletion. Dwell time distributions indicate at least two open and closed states with durations of 0.1 to several seconds, with most channels having an open probability of [&ge;]0.7. We also detect silent channels that colocalize with STIM1 but show no activity over tens of seconds, a population that would be undetectable by whole-cell electrophysiology alone. This method offers an approach to explore CRAC channel gating mechanisms and may be applicable to other Ca2+- permeable channels not amenable to patch-clamp techniques.

BackgroundClinical interpretation of missense variants in the hERG potassium channel encoded by the KCNH2 gene remains a major challenge in inherited arrhythmia syndromes. Functional studies often rely on a minimal set of channel properties, mainly current amplitude measurements, which do not capture the multidimensional nature of channel gating and its impact on ventricular repolarization. We developed a multiscale computational framework to quantitatively link multiparametric channel dysfunction to ECG phenotypes. MethodsWe generated multiparametric electrophysiological profiles for KCNH2 variants using high-throughput patch-clamp, quantifying nine biophysical properties including conductance, voltage dependence, and gating kinetics. These parameters were incorporated into a modified formulation of IKr embedded in a human ventricular electrophysiology model. The resulting framework, termed MeTAL (Multiscale-enriched Transformation and Analysis for Long-QT), produces physiologically calibrated pseudo-ECGs enabling quantitative evaluation of QT dynamics. We systematically analyzed the contribution of individual and combined gating parameters and applied the model to 41 KCNH2 variants across ACMG classes, comparing simulated QTc values with clinical data. ResultsMultiparametric profiling revealed complex functional signatures in most variants, with concurrent gain- and loss-of-function effects affecting distinct gating processes. Simulations demonstrated that ventricular repolarization, though strongly determined by current amplitude, is substantially influenced by inactivation-related parameters, particularly the slope and voltage dependence of inactivation. Interaction analyses showed nonlinear relationships between gating parameters, explaining why variants with similar current density can produce divergent QT phenotypes. In heterozygous simulations, MeTAL reproduced clinically observed QTc distributions across variant classes and accurately predicted the repolarization regime (normal, long-QT, or short-QT) in most cases. ConclusionsMultiparametric integration of ion-channel function within a multiscale electrophysiological model enables mechanistic prediction of QT behavior beyond conductance-based metrics. This approach provides a scalable framework for interpretation of KCNH2 variants and solves the issue of risk stratification in inherited arrhythmia syndromes while offering new opportunities for variant-specific and pharmacological modeling of repolarization.

Apohas Liquid State Intelligence Platform (LSIP) records ellipsometric waveforms from injections depositing sub-microgram quantities of antibody drop-by-drop onto a liquid reservoir. We previously showed that a behavioural feature extracted from the waveforms, VIBE1, identified antibodies carrying multiple biophysical liabilities in an industrial dataset of 71 monoclonal antibodies, and enriched for clinical failure across a larger dataset of 235 therapeutic antibodies [1]. Here, we use an auxiliary coalescence-sensor channel to decode VIBE1 by separating the coalescence event from its propagation through the substrate. The pertitration drop-to-drop standard deviation of pinch-off time,{sigma}{tau} , explains most of VIBE1s variance across the dataset (R2 = 0.92, n = 1182). High-speed imaging at 10,000 frames per second reveals that all imaged drops initially thin at the same Newtonian capillary-inertial rate while the neck remains wide. In drops from certain antibodies, the thinning bridge then decelerates as internal strain builds in the narrowing neck. This elasto-capillary stiffening response has a timescale{lambda} that decreases as pinch-off time{tau} i increases across the imaged set.{sigma}{tau} is therefore a readout of the antibodys propensity to undergo a transient gel-like stiffening response during coalescence, and that variability is what VIBE1 captures. The signal is concentration dependent, and absent in bovine serum albumin (BSA) tested at up to an order of magnitude higher molarity than the antibodies, despite BSA being a strongly surface-active globular protein. The instrument is configured so that complex behaviours of this kind appear in its recorded waveforms; the gel-like coalescence response we identify here is one such phenomenon.

Millimeter-wave radar can quietly monitor health and behavior at home, which is vital for supporting people living with dementia. Most studies, however, remain limited to short-term testing in controlled spaces. Real-world deployment requires robust activity classification as a prerequisite: vital-sign and behavioral sensing require fundamentally different processing pipelines, and absent periods need to be reliably distinguished from stationary states. Bridging the critical gap between controlled laboratory demonstrations and continuous home monitoring, this paper introduces a self-adapting radar framework that extracts meaningful behavioral segments from massive, unconstrained real-world data. The system performs continuous real-time activity classification (stationary, walking, and absent) and target localization, selectively directing downstream processing to the most informative segments. It addresses key real-world deployment challenges including adaptive thresholding across subjects and environments, and walking detection under naturalistic activity conditions. Prior to integration with the Minder platform, the system was validated in a fully instrumented studio apartment against ground truth. Across 12 subjects, the system achieved an overall classification accuracy of 0.98, with F1 scores of 0.99 for absence and stationary states, and 0.95 for walking. Event-based evaluation yielded a per-subject walking sensitivity of 0.916{+/-} 0.058 and F1 score of 0.935 {+/-}0.030. Localization root mean square error during movement was 0.40 m. The results demonstrate reliable performance suitable for transitioning to long-term real-world home deployment.

Developing parsimonious, mechanism-aware quantitative models that predict how biological effectiveness changes with different modifiers remains, in general, an unsolved problem. Advances in radiobiological research have created a large knowledge base of first-principles mechanistic models of radiation response that, in principle, could accurately predict radiosensitivity across different experimental and clinical conditions. However, in practice these mechanistic models come with an overabundance of parameters, the majority of which are practically unidentifiable and, moreover, likely unnecessary if one simply wishes to predict how radiosensitivity changes for some specific modifier of interest. Nevertheless, determining which few details in the full mechanistic model are relevant for a given purpose, as well as how to remove any other extraneous details, remains a highly non-trivial task. In this study, we demonstrate the potential of model reduction, starting from a detailed mechanistic description, as a systematic strategy for deriving parsimonious, experimentally falsifiable radiobiological descriptors. As a proof-of-concept demonstration, we apply the Manifold Boundary Approximation Method (MBAM) to a Mechanistic Model of DNA Repair and Survival (MEDRAS), for the problem of cell survival prediction following an acute exposure. Our findings reveal that the complete MEDRAS model for an arbitrary mixed-quality exposure can be structurally simplified to a reduced three-parameter model for an effective uniform-quality, named MEDRAS-LPL. Additional MBAM analysis on MEDRAS-LPL identifies two boundaries in parameter space, corresponding to sparsely ionizing and densely ionizing radiation. Mapping of MEDRAS-LPL parameter space on to effective LQ space further demonstrates that parameters close to the sparsely ionizing boundary line up with expectations from the theory of dual radiation, while parameters close to the densely ionizing boundary line up with expectations from a purely linear model based on a target-theory description. Moreover, our formalism predicts enhanced synergistic interactions between sparsely ionizing and densely ionizing radiation beyond the Zaider Rossi model (ZRM) paradigm, in line with empirical observations. The results highlight the potential for using reduced-order models not only for predictive applications but also for generating novel hypotheses that can inform future experimental designs and optimization strategies in radiobiology.

This work focuses on the global and partial identification problem for fractional differential equations. We provide a general numerical procedure based on global and local optimization algorithms with two refinements for biological systems that ensure solution positivity and homogeneous parameter units. The method is applied to a new fractional model of Dengue outbreak called the Fractional Homogeneous Nishiura (FHN) model, calibrated using data of newly infected people in Cape Verde. We show that our identification method yields a better fit between data and model solutions than previous approaches and that our FHN model captures the dynamics of Dengue more closely than existing systems.